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ABSTRACT 

The  deformation  of  a  solid  Induced  by  swelling  is  equivalent  to  that  caused 
by  a  temperature  change.  A  generalized  Duhamel- Neumann  form  of  Hooke's 
law  is  employed  to  treat  a  wide  variety  of  environmental  problems  by  the  joint 
application  of  solid  mechanics  and  elementary  physical ohemtstry.  This  approach 
Is  Illustrated  for  a  swollen  fiber  reinforced  material,  employing  physical  chem¬ 
istry  concepts  i  naioromeohanlcs,  and  laminated  anisotropic  plate  theory.  The 
specific  results  are  applicable  to  the  design  of  dimensionally  stable  composite 
materials  in  variable  thermal  or  swelling  environments.  A  new  strain  invariant 
for  laminates  under  these  types  of  environments  is  also  introduced. 

Distribution  of  this  abstract  is  unlimited. 
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INTRODUCTION 


A  general  phenomenon  of  nature  is  the  dilatation  of  a  solid  induced  by  a 
variety  of  causes.  Some  reasons  for  this  phenomenon  are:  an  increase  in 
temperature;  absorption  or  loss  of  swelling  agents  such  as  water,  organic 
solvents,  etc;  and  the  sudden  expansion  of  absorbed  gases  in  a  solid.  These 
natural  processes  are  responsible  for  such  effects  as  thermal  stresses  and 
the  swelling  stresses  (References  1  and  2)  to  be  discussed  in  this  report. 

Swelling  stresses  are  often  involved  in  environmental  cracking  of  solids 
and  the  internal  fracture  of  a  solid  after  it  is  subjected  to  laser  radiation 
focused  on  internal  planes  of  the  material.  In  the  latter  illustration,  stresses 
are  induced  by  the  combined  action  of  locally  large  thermal  expansions  and  the 
rapid  expansion  of  volatile  products  from  the  thermal  decomposition  of  the 
solid.  In  the  everyday  world  these  forces  are  seen  in  the  cracking  of  soil  as 
it  dries,  the  expansion  of  door  jambs  and  window  sashes  in  humid  weather,  and 
the  cracking  of  paint  as  it  dries.  Figure  1  is  a  photograph  of  a  rubber  block 
which  was  swollen  and  subsequently  dried  out,  primarily  by  nonuniform 
evaporation.  The  resulting  stresses  fractured  the  block  in  a  cracking  pattern 
similar  to  that  observed  in  dried  soil  or  on  a  painted  surface. 

Despite  the  considerable  work  devoted  to  thermal  expansion  and  thermal 
stress  problems,  the  current  literature  does  not  provide  a  general  basis  for 
formulating  an  analysis  of  such  diverse  yet  related  phenomena  as  outlined 
above.  The  Importance  to  modem  engineering  technology  and  the  biological 
sciences,  however,  has  prompted  this  attempt  at  an  analysis  of  the  conse¬ 
quences  of  such  phenomena.  For  the  sake  of  brevity  and  laboratory  conve¬ 
nience,  we  shall  analyze  the  phenomenon  of  swelling  as  an  illustrative  example 
and  show  how  this  problem  parallels  the  related  thermal  problem. 
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EXPANSIONAL  STRAINS 


In  this  section  we  shall  establish  the  procedures  for  the  computation  of 
strain  when  an  object  Is  deformed  by  swelling  or  Internal  gas  expansion.  With 
swelling  we  are  concerned  with  the  equilibrium  existing  between  at  least  two 
phases.  The  simplest  case  to  consider  is  that  in  which  one  phase  contains  two 
components,  while  the  other  contains  only  one.  We  shall  have  a  particular  inter¬ 
est  in  the  specific  case  where  the  mixed  phase  is  solid  and  the  pure  phase 
liquid. 

The  absorption  of  a  vapor  or  liquid  by  a  solid  is  a  natural  process,  not 
uniquely  dependent  on  any  specific  attraction  between  the  solid  in  question  and 
the  swelling  agent,  but  rather  on  the  general  diffusing  tendency  of  two  sets  of 
molecules.  This  diffusing  tendency  is  no  different,  in  principle,  from  that 
existing  between  two  liquids  of  similar  ohemioal  constitution:  for  example,  the 
mixing  of  alcohol  and  water.  These  remarks  are  particularly  true  for  organic 
polymeric  materials  absorbing  organic  vapors  and  fluids,  whioh  are  predomi¬ 
nantly  entropy  controlled  processes.  However,  for  some  water-swelling  mate¬ 
rials,  the  driving  force  in  the  swelling  prooess  is  the  chemical  attraction  be¬ 
tween  water  and  solid,  which  is  effective  inBpite  of  an  adverse  entropy  reduction. 
Detailed  discussions  on  these  points  are  available  in  standard  tests  (References 
3  -  0)  on  theory  of  solutions. 

The  extent  of  swelling  will  be  defined  by  the  volume  fraotion  v2  of  solid 

in  the  mixture  of  solid  and  liquid.  On  absorbing  a  quantity  of  vapor  or  liquid  a 

unit  isotropic  cube  in  the  dry  state  deforms  into  a  cube  of  edge  X.  The  volume 

3 

swelling  ratio  for  an  incompressible  solid  will  then  be  given  by  l/v2  or  \  . 

The  normal  strain  corresponding  to  this  volume  change  is  defined  as 

€  =  A  “I  s  (  I  /v2  -I  (I) 

for  small  strains,  <  .  These  definitions  are  independent  of  the  nature  of  the 
swelling  agent  or  of  the  question  of  equilibrium  between  the  specific  fluid  and 
swollen  solid.  We  consider  v2  as  a  parameter  which  defines  the  state  of  strain 
in  an  unstressed  swollen  Infinitesimal  volume  element.  It  is  sufficient  at  this 
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point  to  simply  recognize  that  procedures  exist  In  the  physical  chemistry  of 
mixtures  and  solutions  (Reference  3)  which  allow  for  the  nreHirHnn  nf for 
specific  combinations  of  fluids  and  solids. 


During  the  performance  of  our  conceptual  experiment  tho  isotropic  body 
was  observed  to  undergo  a  pure  volume  dilatation  in  the  absence  of  surface 
tractions.  If  it  were  only  required  that  we  pass  from  the  initial  volume  to  the 
final  volume  under  no  surface  tractions,  the  same  state  of  hydrostatic  strain, 
<,  could  have  been  produced  by  a  simple  temperature  rise.  A  thermal  expan¬ 
sion  is  therefore  the  mechanical  equivalent  of  a  swelling  process. 


We  are  now  in  a  position  to  consider  the  form  of  the  constitutive  relations 
In  a  more  general  environmental  state.  Performing  the  well  known  power 
series  expansion  of  the  strain  energy  function  we  have 


ZW  =  Cc  +Cjj  e..  +C,jkl  e,j  ek|  +  — - 


(2) 


where  the  coefficient-’  CQ,  CJj,  etc.  are  functions  of  the  entropy,  the  tempera¬ 
ture,  and  such  composition  variables  as  v0.  Since  W  *  0  when  =  0  by  defini¬ 
tion,  the  constant  C  must  vanish.  Furthermore,  since 


°ii  T  cii  +  0<V  151 

Cjj  must  be  equal  to  20^  at  zero  strain.  Such  stresses  occur  in  a  solid  that  is 
completely  constrained  against  deformations  caused  by  environmental  factors 
such  as  swelling  and  change  in  temperature.  In  what  follows  we  shall  call  these 
deformations  expansional  strains  The  simplest  assumption  to  treat  these 
effects  is  to  allow  Cjj  to  be  proportional  to  the  expansional  strums,  i.e. , 


Clj  s-2Cijkl  * kl  (4) 

where  s  a^j  AT  for  temperature  changes  and  to  employ  Equation  1  for  cal¬ 
culating  expansional  strains  due  to  swelling  or  internal  gas  expansion.  When 
we  drop  higher  than  second  order  terms  in  the  strain  energy  density,  we  get 

W  s"Cijkl  *lj  €ki  +  \  Cijki  eij  ®ki  (5) 
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Hence,  the  stress  components  are  given  by 

c r,.  *  C...  (e.  )  (6) 

ij  ijkl  kl  kl 

which  is  a  more  genera)  version  of  the  Duhamel-Neumann  form  of  Hooke's  Law. 
The  inverse  of  Equation  6  is  then 

* i|  1  Sijkl  ^kl  +  *ij  {71 

where  Sjjy  is  the  compliance  tensor. 

Since  the  expanslonal  strains  are  the  components  of  a  second  order  tensor, 
they  must  obey  the  tensor  transformation  law 

*  ij  "ki  'u  *w 

where  tjj  are  the  direction  cosines  between  xj  and  xj  .  For  the  specific  case  of 
an  orthotropic  material,  where  unprimed  coordinates  refer  to  the  axes  of  ma¬ 
terial  symmetry,  the  expanslonal  shear  strains  are  zero  and  Equation  8  reduces 
to 

*i|*'n  'n  *n  *«, *„  (9) 

An  experimental  verification  of  this  transformation  law  will  be  given  in  a  sub¬ 
sequent  section. 


For  an  isotropic  solid,  Equation  7  reduces  to 

*lj  *  2G~  °ij  66  ”  9iT,^lj  +<^(j 

where  8(j  is  the  Kronecker  delta,  G  the  shear  modulus,  K  the  bulk  modulus, 9 
=  <x.  ,  and  «  is  the  linear  expanslonal  strain  defined  by  Equation  1  for  swelling 
and  by  qAT  for  thermal  expansion.  In  the  general  case,  small  expanslonal 
strains  from  different  sources  are  additive  so  that 

*  €Thermol  +<Swellinq  + 

The  corresponding  strain  energy  function  is  given  by 


I ' 
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Accordingly,  we  mnv  ennelttHe  this  wMinn  h„  n^ti  M»<»f  n  n*l  MViinl*nt  <nw» 

-  -  -  -  •  *  - -  "«/  - - *0  “w*  ***v  J/»*JW*V**»  p*  V 

cess  which  yields  an  expansional  strain  can  be  formulated  into  a  boundary  value 
problem  within  the  framework  of  formal  elasticity  theory.  In  fact,  solutions 
obtained  in  thermal  elasticity  can  be  extended  to  general  expansional  problems 
by  the  substitution  of  «tj  for  aij  AT.  This  statement  should  also  be  valid  not  only 
for  uniform  expansional  fields  but  also  for  problems  where  expansional  gradients 
are  of  concern.  Where  tbo  properties  of  a  thermal  gradient  and  thermal  diffu¬ 
sion  are  expressed  in  Fourier's  law  (Reference  7),  the  problem  of  gas,  vapor, 
or  liquid  diffusion  and  concentration  gradient  will  be  governed  by  Flok's  law 
(References  1  and  8). 

A  problem  in  expansional  elasticity  shall  be  discussed  and  tested  experi¬ 
mentally  in  the  next  section.  The  problem  is  to  establish  as  a  function  of  con¬ 
stituent  material  properties  the  expansional  strains  in  a  thin  heterogeneous 
sheet  submerged  in  a  swelling  liquid  and  then  to  predict  the  swelling  properties 
of  an  angle-ply  laminate  constructed  of  this  heterogeneous  material. 

The  problem  of  an  advancing  boundary  of  swollen  material  within  a  homo¬ 
geneous  isotropic  cylindrical  body  has  been  treated  by  Alfrey,  Gurnee,  and 
Lloyd  (Reference  1)  by  use  of  elasticity  theory.  Our  prime  concern  in  subse¬ 
quent  sections,  however,  will  be  to  study  the  influence  of  swelling  in  hetero¬ 
genous  anisotropic  systems. 
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ANALYSES  Or  A  UNIDIRECTIONAL  COMPOSITE  SHEET 


The  prediction  of  the  elastic  and  thermal  properties  of  unidirectional  com¬ 
posites  as  functions  of  constituent  material  properties  and  phase  geometry  has 
received  considerable  attention  in  recent  years  (References  9  and  10).  Very 
recently  Sohapery  (Reference  11)  has  derived  upper  and  lower  bounds  aa  well 
as  convenient  approximate  expressions  for  thermal  expansion  coefficients  based 
on  the  principles  of  minimum  potential  and  complementary  energy.  We  shall 
now  demonstrate  that  Schapery's  solutions  are  equally  valid  for  expansional 
strains  induced  by  a  swelling  agent. 

Because  our  argument  is  dependent  on  experimental  verification  wa  have 
chosen,  for  experimental  ease,  an  elastomeric  matrix  material  reinforced 
with  a  system  of  parallel  nylon  textile  cords  for  the  demonstration.  The  me¬ 
chanical  characterisation  and  the  special  experimental  techniques  required  for 
orthotropic  materials  has  previously  been  reported  in  References  12  and  13. 

The  constituent  materials  properties  are  as  given  below,  where  the  sub¬ 
script  f  stands  for  fiber  and  m  for  matrix: 

Ef  *  292,000  psi 
E_  *  300  psi 

ffl 

vf  *  0.2 
v_  *  0.4999 

fn 

The  composite  moduli  below  were  determined  experimentally  and  are  in  agree¬ 
ment  with  the  ndcromechanicB  prediction  of  Hermans  (Reference  14): 

E„  .  I  32,000  psi 

Et2  *  1050  psi 

p(t  *  0.36 

G,c  >  263  psi 

where  Ell  and  E22  are  the  respective  Young's  moduli  parallel  and  normal  to 
the  fiber  direction,  v12  is  Poisson's  ratio  measuring  transverse  contraction 
under  a  uniaxial  stress  acting  parallel  to  the  fibers,  and  G12  is  lbs  shear 
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modulus  with  respect  to  the  axes  parallel  and  normal  to  the  fibers.  Owing  to  its 
high  degree  of  anisotropy,  this  material  offers  a  stern  test  for  the  hypothesis 
presented  above. 


We  shall  first  offer  a  graphical  demonstration  of  the  transformation  prop¬ 
erty  of  C|j ,  l.e. ,  Equation  9.  This  was  achieved  by  immersing  unidirectional 
sheets  possessing  different  fiber  orientations  In  an  effiolent  swelling  agent  for 
the  matrix,  namely  benzene.  Results  of  a  typical  experiment  are  shown  in 
Figure  2(a)  for  a  30-degree  specimen  in  the  dry  and  swollen  states. 


In  the  performance  of  the  experiment  the  specimens  were  observed  to 
undergo  rather  large  strains  invalidating  the  assumption  of  infinitesimal  de¬ 
formations.  Accordingly,  to  obtain  a  precise  comparison  of  theory  and  experi¬ 
ment,  one  must  employ  a  geometrically  nonlinear  measure  of  strain.  We  shall 
utilize  the  Lagranglan  or  Green  strain  tensor  defined  (Reference  7)  as 


with 


2« 


jk 


1  V  +  “M  * 


“l.i  “l.k 


The  axtenaional  strains  were  computed  as 

[(E.  +  i>8-'] 


and 


[lE.+n* 


(13) 


(14) 


where  E  *  Al/1. 

< 


The  corresponding  shear  strains  were  computed  as 
•*  *  -f[  0  +  2«i,  >  (■ 


(15) 


whBie  thB  change  of  the  angle  between  two  originally  orthogonal  coordinate  axes 
is  given  as  K  was  then  assumed  that  a  linear  relationship  exists  between 
the  stresses  (baaed  upon  undeformed  dimensions)  and  the  Lagranglan  strain 
components.  TUb  assumption  is  Justified  only  In  that  it  leads  to  a  successful 


-v- :  .-aUimiAlil-u.-.  .  iUr.aiwA:  yj  H  u- — • 
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Figure  2.  Comjiarison  of  Dry  and  Swollen  Anisotropic  Specimens: 

a)  30-degree  Unidirectional  Materiali  and 

b)  An  Angle-Ply  Laminate 
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correlation  between  the  analytlcaKknd  experimental  results  for  the  specific 
materials  an-i  problem  considered  here.  The  choice  of  a  strain  tensor  depends 
upon  the  magnitude  of  the  experimental  effects  and  is  incidental  to  the  theme  of 
this  presentation. 


In  the  plane  of  a  thin  sheet.  Equation  9  reduces  to 


+  n8  « 

II 

II 

ta 

c 

=  n  € 

+  m8 

«* 

as 

€lt 

«  mn  ( -  «  „  ) 

where  the  xj  axis  is  parallel  to  the  fibers  and  is  parallel  to  the  long  edge  of 
the  sheet*  m  «  oos  9  and  n  »  sin  8.  To  illustrate  the  transformation  of  Equations 
16,  the  expansional  strains  «Jj  were  measured  and  computed  in  accordance 
with  Equations  14  and  IS  as  a  function  of  8,  the  data  for  being  indicated  in 
Figure  3.  The  solid  eurve  signified  by  the  term  unidirectional  represents  the 
first  transformation  equation  of  (18)  with  s  0.02  and  «22  *  °-75-  The  shear 
component  e|2  ,B  presented  in  Figure  4  and  ia  in  accord  with  the  above  cited 
transformation  rules. 


i 

f 

\ 


To  attempt  the  prediction  of  the  composite  expansional  strains  based  upon 
Its  constituent  material  parameters,  the  expansional  strains  of  the  fiber  and 
matrix,  «f  and  «m  respectively,  In  benzene  were  determined  by  experiments  as 

•_  «  0.84 

m 

«f  «o.oi 

These  values  as  well  as  the  moduli  presented  earlier  were  used  to  compute 
the  composite  expansional  strains  and  «22  from  Sohapery's  formulas  which 
were  developed  for  linear  thermal  expansion  coefficients.  For  the  longitudinal 
expansional  strain  we  have 


« 


ii 


,  Em  vm  +  Ef«fwf 
Em  vm  +  Efwt 


1 17) 


hh. '  XaV . :  ■  -.‘“jtv  %V  **Ak«'  fl iijiA \U» *•'  «»V>  - iuJ •*  jJBfcia  '  ' 
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NQLE  PLY 


EXPANSIONAL  SHEAR  STRAIN 
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Figure  4.  Comparison  of  Theoretical  and  Experimentally  Observed 
Expansional  Shear  Strains  for  a  Unidirectional  Material 
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where  v  stands  for  volume  fraction.  For  a  volume  fraction  of  fibers  equal  to 
0.45  the  use  of  Equation  17  gives 

«„  =0.011 

The  expansional  strain  normal  to  the  fiber  direction  is  given  as 

**2  s(  l  +  vm  >  +  <l  +  V1  )€fvf  <vfvf*  vmvm)  <l8> 

which,  In  our  composite,  yields  a  valu  '  of  0. 70.  Both  of  these  estimates  are 
well  within  the  limits  of  our  experimental  precision.  The  swelling  expansion 
for  this  matrix  material,  «m»  Is  In  turn  predictable  from  the  well  known 
Flory- Huggins  equation  In  the  polymer  physical  chemistry  literature  (Refer¬ 
ences  3,  5,  and  6).  Thus  we  are  able  to  form  a  continuous  development  for  the 
calculation  of  a  physical  property  ranging  from  file  physical  chemistry  of 
constituent  materials  to  the  mechanical  response  of  a  composite  material.  In 
the  next  section  we  shall  extend  the  theoretical  chain  by  predicting  the  expan- 
sional  properties  of  a  laminate  from  the  knowledge  of  the  properties  of  the 
individual  layers. 
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Onoe  the  expans lonal  properties  of  a  single  layer  of  fiber  reinforced 
material  have  been  computed  from  the  constituent  material  properties  as 
discussed  in  the  previous  section,  the  response  of  a  laminate  (a  system  of 
unidirectional  layers  of  various  orientations)  can  be  determined  through  the 
use  of  classical  laminated  plate  theory  (Reference  15).  In  this  theory,  each 
layer  Is  treated  as  a  homogeneous  anisotropic  material  In  a  state  of  plane 
stress.  Hence  the  constitutive  equations  for  an  individual  layer  are  of  the  form 

ai  a°i}  <#i  ~ V*  l,j  “  *‘2’6  <■») 

where  standard  contracted  notation  (Reference  15)  is  now  being  employed* 
l.e.  *  Qjj  representing  the  reduced  stiffness  matrix  for  plane  stress.  ©J  the 
engineering  strain  components,  and«j  the  expanalonal  strains.  In  Equation  19 
the  normal  components  of  stress  and  strain  are  represented  by  the  subscripts 
1  and  2*  while  the  shear  components  are  indicated  by  the  subscript  6.  The 
stress  resultants  Nj  are  then  given  by 

»  h/t 

N,  *  J  0{j  <Sj  -  «j)  d«  (20) 


where  h  Is  the  plate  thickness.  For  laminates  in  which  no  bending  effects  are 
present*  the  strains  6j  are  constants;  hence  when  Nj  -  0,  we  get 


A«  *j  * 


where 


Vj  " 


a,,  « 


Equation  21  defines  the  laminate  strains  under  pure  swelling  or  thermal 
environments*  and  therefore  the  equivalent  laminate  expansional  strains  c  j, 
<2*  and  «g.  In  what  follows*  we  shall  restrict  our  attention  to  uncoupled  lam¬ 
inates  consisting  of  layers  of  the  same  material  and  for  which  Alg  =  A2g  =  0, 
so  that «  g  also  vanishes.  An  example  of  such  a  laminate  is  a  plate  composed 
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of  a  series  on  angle-ply  systems  stacked  symmetrically  with  respect  to  the 
central  plane.  Using  the  invariant  properties  of  the  reduced  stiffness  transfor¬ 
mation  equations  (Reference  15),  the  expans  tonal  strains  can  be  oonviently 
expressed  as 

„  .  Ag2  R|  ~~  A|g  ^2 

•  ~  a  a  rAi 

(i  Mta  Mia 

(2  3] 

e  s  AlL_R.1.1-AJj,ll- 

8  A  A  —  A* 

Mu  Mia 

where 

R,  «  J,  h  +  Ja  H, 

Rg  a  J|  h  —  Jj  H| 

A,,  «  U,  h  +Ua  H,  +U3  H2 
Aaa  *  U,  h  -  u2  H,  +  U,  H2 
A  n  ■  U4h  “  U8  Hg 
J,  •  (U,  +  U4  )  W.  +  2U.  W. 

( 24 

jf  .  uaw,  +2w,  ( u |  +eu,-u4 ) 

«,  ■  •£<«]-  +««  ) 

w«  -  -i-  > 

H,  «  I  hnqos2  0n 
ml 
N 

Ht  *  t  hn  cos  4 Bn 
ml 

and  and  «  2  are  the  layer  expanslonal  strains  parallel  and  normal  to  the 

fibers,  respectively,  h  Is  the  laminate  thickness,  N  the  number  of  layers, 

_  th 

an  the  angle  between  the  x1  axis  and  the  fibers  in  the  n  layer,  and  h„  the 
til 

thickness  of  the  n  layer.  The  various  Uj  ore  functions  of  the  orthotroplo 
elastic  moduli  of  the  layers  and  are  defined  In  Reference  15. 

Previously  reported  material  properties  were  used  to  compute  the  lami¬ 
nate  expanslonal  strains  for  angle  ply  composites  and  are  shown  as  the  "angle 
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ply1'  curve  in  Figure  3.  These  calculations  may  be  somewhat  surprising  in  that 
iiasy  suggest  that  at  certain  angles  •  the  laminate  will  contract  in  the  longitudinal 
direction  when  heated  or  swollen  Figure  2.  The  triangular  data  points  were  ob¬ 
tained  from  equilibrium  swelling  of  various  angle-ply  samples.  As  indicated  in 
the  figure  we  have  observed  a  10%  decrease  in  length  for  a  ±15-degree  laminate. 
Note  that  the  15.de  gree  unidirectional  material  gave  a  corresponding  positive 
expansion  of  15%.  Dow  and  Rosen  (Reference  16)  are  observing  this  same  effect 
in  thermal  expansion  problems.  The  obvious  agreement  between  theory  and 
experiment  lends  confidence  in  the  approach  outlined  here.  Calculations  of  this 
type  will  have  technological  significance  as  they  specify  the  material  parameters 
and  geometric  construction  necessary  to  yield  structural  elements  with  zero  or 
minimal  expansions!  coefficients. 
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ISOTROPIC  EXPANSIONAL  STRAIN  IN  A  LAMINATE 


Consider  the  problem  of  defining  a  lamination  sequence  such  that  the  in¬ 
plane  expanslonal  strain  field  is  isotropic  •  i.  e. ,  the  expansional  normal  Btrain 
is  the  same  in  every  direction  in  the  plane  of  the  sheet,  while  the  shear  strain 
vanishes.  A  necessary  and  sufficient  condition  for  this  situation  is  that  *  «2. 
In  Equation  23  this  occurs  when  the  following  relation  is  satisfied 

Ri  Ah  4 An 

In  general,  the  isotropic  strain  induced  when  Equation  26  is  invoked  is  a  function 
of  the  details  of  the  lamination  sequence.  However,  if  we  consider  :he  case 
where  A^  *  A22,  we  find  from  Equations  24  and  25  that  F  0,  which  implies 
that  Rj  *  R2  1  J  jh.  If  this  is  true  we  can  use  Equations  •  and  24  to  obtain  the 
following  expression  for  the  isotropic  strain 

«g«.  s  «gs  -U~,  4~U4  "  (Z6) 

which  is  an  invariant  quantity.  The  expression  in  Equation  26  depends  solely  on 
the  mechanical  and  expansional  properties  of  a  unidirectional  sheet  of  the  lami¬ 
nate  material.  We  can  therefore  assert  the  following  general  principle:  For  any 
uncoupled  laminate  with  equal  stiffnesses  in  two  in-plane  directions  (Aj^  :  A??), 
the  expansional  strain  field  is  isotropic  in  this  plane  and  the  normal  strain  in  all 
directions  is  simply  given  by  Equation  26.  This  general  statement  is  based  on 
the  premise  that  all  layers  consist  of  the  same  material. 

The  statement  holds  despite  the  fact  that  the  laminate  is  not  quasi-isotropic 
with  respect  to  ln-plane  stiffness.  Examples  of  laminates  which  undergo  iso¬ 
tropic  expansional  strain  are  0  -  90*degr«*e  bidirectional  composites,  ±4 5 -degree 
angle  ply,  and  the  combined  angle  ply+a,  -cl,  (ir/2  +0 ),  (•jr/2-a)  for  any  value 
of  a,  all  of  which  consist  of  layers  of  equal  thickness.  For  a  given  composite 
material,  all  of  these  systems  undergo  identical  expansional  strains.  The  par¬ 
ticular  case  of  a  20,  -20,  70,  110-degree  laminate  was  investigated  experi¬ 
mentally  (Figure  5).  In  this  experiment  an  isotropic  expansional  strain  of  0.  034 
was  observed;  the  prediction  was  fj  :  «2  s  0-027. 
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(b) 

Figure  5.  Comparison  of  a)  Dry  and  b)  Swollen  Composite  Materials: 
20,  -20,  70,  ancKlO  Degree  Laminates 
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Equation  26  may  also  be  expressed  in  terms  of  engineering  constants  as 

2  i  E„  -Ejt  )#a 

«  *  W.  +  -  — r - TT -  (27 

1  E»i  +  <*  +  2*rt>  E« 

where  all  moduli  refer  to  the  properties  of  the  unidirectional  material  and  Wp 
Wj  are  given  in  Equations 24.  According  to  Equation  27,  the  isotropic  strain 
becomes 


«  *  W,  +  2W, 


«,  ,09 


Thus,  for  very  highly  anisotropic  composites  l.e.,  high  ratio  of  Eii/Egg)*  ^ 
isotropic  strain  approaches  the  expansional  strain  of  a  unidirectional  ply  paral¬ 
lel  to  the  filaments.  This  same  limiting  oase  may  be  approached  by  employing 
ribbon  shaped  reinforcements  (Reference  19)  possessing  high  aspect  ratios. 
Accordingly,  the  material  fabricator  has  several  means  at  Us  disposal  for 
creating  a  dimensionally  stable  material  for  technological  applications. 


Under  pure  thermal  or  swelling  environments,  the  boundaries  of  a  laminate 
are  stress-free;  however,  each  layer  is  subjected  to  stresses  which  can  become 
quite  large.  In  fact,  in  our  swelling  experiments,  the  induced  stresses  in  the 
laminates  were  frequently  large  enough  to  cause  delamination  and  fracture  of 
the  material.  These  internal-ply  stresses  can  be  estimated  from  the  equations 
of  classical  laminated  plate  theory, 


where 


«<n)  ; 
a 


-cjn>)  for  i,j  •  1,2,6 


+  2  Wt  cos  2  0n 
-  2  cos  2 9n 


-  4W, 


►  in  2  0. 


Angle-ply  thermal  expansion  coefficients  for  various  composite  materials 
are  shown  In  Figure  6. 
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Figure  6.  ComimrUon  of  Hworettoel  Linear  Itarmal  Expanelon 
Coe  (Hole  Ota  for  Several  Structural  Conpoettee 
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CONCLUSIONS 

We  have  shown  here  that  the  changes  in  sice  and  shape  induced  by  the 
swelling  of  a  solid  produces  effects  equivalent  to  those  caused  by  temperature 
changes.  The  generalized  Duhamel- Neumann  form  of  Hooke's  Law  can  be  em¬ 
ployed  to  treat  a  wide  variety  of  environmental  problems  by  the  Joint  application 
of  solid  mechanics  and  elementary  physical  chemistry.  While  we  have  limited 
ourselves  to  equilibrium  problems  ,  these  results  may  be  extended  to  time 
dependent  phenomena  through  the  structure  of  linear  viscoelasticity  (References 
17  and  18)  and  the  proper  definition  (Reference  17)  of  a  new  reduced  time  para¬ 
meter  oc  to  account  for  the  effects  of  the  swelling  agent  on  die  time  scale.  The 
reduced  time,  £ ,  for  combined  thermal  expansion  and  swelling  will  be  defined 
as 


as  discussed  In  Reference  17. 
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